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This note presents neutron scattering cross sections for polyethylene at 296 K, 77 K and 4 K derived from a 
new scattering kernel for neutron scattering off of hydrogen in polyethylene. The kernel was developed in ENDF- 
6 format as a set of S{a,f3) tables using the LEAPR module of the NJOY94 code package. The polyethylene 
density of states (from to sub eV) adopted to derive the new kernel is presented. We compare our calculated 
room temperature total scattering cross sections and double differential cross sections at 232 meV at various 
angles with the available experimental data (at room temperature), and then extrapolate the calculations to 
lower temperatures (77K and 4K) . The new temperature dependent scattering kernel gives a good quantitative fit 
to the available room temperature data and has a temperature dependence that is qualitatively consistent with 
thermodynamics . 



1. INTRODUCTION 

A prototype solid deuterium (SD2) Ultra Cold 
Neutron (UCN) source has been under develop- 
ment for testing at the Los Alamos Neutron Sci- 
ence Center (LANSCE) since the summer of 1998. 
This involves the use of a cold "moderator" sur- 
rounded by a beryllium flux-trap in which is em- 
bedded a tungsten spallation target. The pur- 
pose of the "moderator" /flux-trap system is to 
optimize a flux of cold (20 to 40 K) neutrons 
that interact with a volume of SD2 at 4 K to 
produce UCN through inelastic scattering. Dur- 
ing the prototype testing phase, in order to avoid 
safety issues associated with the use of liquid deu- 
terium or hydrogen, we have chosen polyethylene 
at 4 K as the cold "moderator" material. 

In order to simulate the neutronic behavior of 
this prototype with Monte Carlo methods using 
MCNP, it was necessary to develop a scattering 
kernel in ENDF-6 format for incoherent neu- 
tron scattering from hydrogen in polyethylene at 
cryogenic temperatures, which are not yet avail- 
able in the official neutron scattering database. 
This report describes the development of such a 
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kernel in the form of S(a, /3) tables for polyethy- 
lene at 296 K, 77 K, and 4 K, temperatures rele- 
vant to the different cooling stages of the experi- 
ment. 



2. POLYETHYLENE DENSITY OF 
STATES 

A scattering kernel for neutrons scattering off 
of the hydrogen in room temperature polyethy- 
lene was developed a number of years ago by 
Sprevak and Koppel|21 using the GASKETPm 
code and applying the incoherent approximation 
.The incoherent approximation is appropriate for 
polyethylene since the dominant scattering cen- 
ter is the H nucleus where the coherent cross sec- 
tion accounts for only 2.2% of the total scattering 
cross section {atot — 81.67 barn). The density of 
states of polyethylene used for this scattering ker- 
nel calculation was derived from the work of Lin 
and Koenig[S], who calculated the structure force 
constants assuming polyethylene to be an infinite 
chain of CH2 radicals (methylene) . The force ma- 
trix is thus simplified into finite dimensions as a 
single methylene plus a relative phase parameter 
between adjacent methylene groups. The density 
of states (referred to as the frequency distribution 
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Figure 1. Density of states for polyethylene used 
for the scattering kernel. 
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Figure 2. Sub-thermal region of the polyethylene 
density of states. 



in some literature) defined as 

can be calculated once the energy spectrum tOg is 
determined. Here s indexes energy modes. The 
intra-molecular excitations include bond stretch- 
ing (symmetric and antisymmetric modes), wag- 
ging, bending and rocking 6 . They all contribute 
to the optical modes in the density of states shown 
in Figure ^ Direct excitations of bond stretch- 
ing results in a sharp peak at 0.36 eV, and the 
other relative motions between bonds contribute 
to excitations between 0.08 and 0.19 eV. 

The sub-thermal energy part of the spectrum 
results from the collective excitations of the solid 
lattice, which can be approximated by the Debye 
spectrumj2j for phonon excitations. A number 
of years later, Swaminathan and TewarijT] de- 
veloped a more detailed polyethylene density of 
states in the sub-thermal energy range as part of 
their efforts to explore the use of cold polyethy- 
lene as a source of cold neutrons and to study 
the effect of the degree of crystallinity on neutron 
scattering. They proposed an anisotropic disper- 
sive continuum model for the acoustic phonons 
in crystalline polyethylene and a density of states 
as shown in Figure |21 Their paper [7] contains a 
brief description of some of the other scattering 
kernels that have been developed for polyethylene 
over the years. 



Our UCN source involves a broad spectrum of 
neutrons, from the spallation energy of several 
tens of MeV to a few hundred neV. To perform 
accurate simulations we require a good knowl- 
edge of our neutron moderator, polyethylene, over 
most of this energy range. We approach this 
by combining the epithermal part of the Sprevak 
and KoppeiPl spectrum with the low energy part 
of the Swaminathan and Tewari^ spectrum for 
crystalline polyethylene. The resulting polyethy- 
lene spectrum is shown in Figure Q The low en- 
ergy (< 30 meV) component of the spectrum, due 
solely to the acoustic phonons, should be normal- 
ized to properly account for the number of de- 
grees of freedom (3 translational for each of the 
N scatterers(CH2)) in a volume V; i.e., 

/ dujp{uj) = iN/V. (2) 

Here Tq is the Debye temperature of the 
polyethylene. This combined spectrum was 
then run through the LEAPR 8| module of the 
code package to obtain scattering ker- 
nels at various temperatures in ENDF-6 format. 

3. LEAPR FORMALISM 

Temperature dependent inelastic scattering 
components can be calculated by integrating the 
double-differential cross sections obtained from 
the scattering kernel; i.e. the S(a,/9) tables, over 
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energy and angle: 

(fa{T,E, -> Ef,e 
dEfdn 



^(«,/3),(3) 



where ab is the bound scattering cross section. 
The unit-less momentum transfer a and energy 
transfer /3 are defined as 



Ei + Ef - 2 y/E^Ef cost 
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/3 
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where A is the atomic weight. The scat- 
tering kernel S{a^(3) for inelastic phonon cre- 
ation/annihilation is 



1 1"^ 



(4) 



with 
7(i) = a 



2/3sinh(/3/2) 



(1 



)e->^/^df3 (5) 



incorporating the material characteristics 
through its density of states p{f3)). 

The incoherent elastic component has a simpler 
formalism: 

I _ -iEW X 

-iE)=<r,^^ and W=—. (6) 

where the Debye- Waller factor W is determined 
once the Debye- Waller coefficient A is calculated 
by LEAPR. 

4. COMPARISON OF THE SCATTER- 
ING KERNEL WITH EXPERIMEN- 
TAL DATA. 

We have tested the proposed polyethylene den- 
sity of states by using it to generate room tem- 
perature cross sections in LEAPR and comparing 
the results with available data. A comparison of 
the calculated and measured [Til] total scattering 
cross section for room temperature polyethylene 
is shown in Figure|31 The experimental data have 
been adjusted by subtracting a 1/v hydrogen ab- 
sorption cross section (equal to 0.334 b/H atom 
at 0.0253 eV) and a constant 4.74 barns for scat- 
tering from the carbon nucleus. The agreement 
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Figure 3. Total neutron scattering cross section 
from room temperature polyethylene, calculated 
using the density of states presented in FigQ] 
The result of calculation (solid curve) is compared 
with experimental data (dots). 



between the theory and the measurements is ex- 
cellent. 

The polyethylene density of states is further 
tested by looking into the details of the double- 
differential cross-sections obtained from LEAPR. 
A direct comparison of these double-differential 
cross sections with experimental data^J at three 
different angles for an incident neutron energy of 
232 meV (T„ = 2700 K) is given in Figure H 

The experimental data were taken from figures 
in the Sprevak and Koppel paperl^. Once again, 
the agreement is satisfactory. Overall, the struc- 
ture built into the S(q;,/?) tables from the pie) 
matches well the structure observed in the neu- 
tron scattering. The peaks at low energy cor- 
respond to the excitations of the energy states 
around 0.16 eV and 0.095 eV (seen in Figure 
In addition, there is little up-scattering (no 
energy-gain peaks). This is what we would expect 
in situations of predominately inelastic scattering 
where the incident neutron energy is much larger 
than the moderator energy (T„ >> T). 
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Figure 4. Double differential cross sections for 
232 meV incident neutron at 25, 60 and 90 de- 
gree scattering off of room temperature (296K) 
polyethylene. The results of calculation (solid 
curves) are compared with experimental data 
(dots). 
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Figure 5. Calculated neutron elastic scattering 
cross section for polyethylene temperatures of 
296, 77 and 4 K. Also shown are the Debye- Waller 
coefficients at each temperature. 



5. EFFECT OF POLYETHYLENE TEM- 
PERATURE ON COLD NEUTRON 
SCATTERING. 

Equation |^ implies that the elastic scattering 
cross section is temperature independent in the 
limit of low neutron energy and inversely depen- 
dent on temperature in the high energy limit. 
Plotted in Figure O is the effect of polyethylene 
temperature on incoherent elastic scattering. As 
expected, the elastic scattering of cold neutrons 
(T„ < lOK, or 1 meV) is independent of polyethy- 
lene temperature, and the cross section at high 
energies scales roughly as 1/T. 

For incoherent inelastic scattering, on the other 
hand, Figure |^ shows that low energy neutron 
scattering is strongly dependent on polyethy- 
lene temperatures while epithermal scattering 
(T„ >1000 K or 100 meV) is almost tempera- 
ture independent. In the high energy region (rela- 
tive to the moderator temperature) the dominant 
process is down-scattering through excitation of 
internal states of the material. These states are 
intrinsic and their availability depends very little 
on temperature. For neutrons at very low ener- 
gies, on the other hand, the thermalization pro- 
cess requires up-scattering where the cross sec- 




Figure 6. Calculated neutron inelastic scatter- 
ing cross section for polyethylene temperatures 
of 296, 77 and 4 K, compared with the nuclear 
absorption cross section. 



tion scales as the number of available heat carri- 
ers; i.e., the phonon density. The phonon number 
density is governed by the Boltzmann factor and 
can be suppressed by cooling the material. These 
solid state effects produce the observed tempera- 
ture dependence. 

Moreover, it should be noted that the useful- 
ness of CH2 at 4 K as a moderator is limited by 
nuclear absorption. The absorption cross section 
scales inversely with the neutron velocity. For 
CH2 at 40 K, inelastic scattering has a larger am- 
plitude than absorption at all neutron energies. 
For CH2 at 4 K, on the other hand, the nuclear 
absorption overwhelms the inelastic scattering for 
all neutrons with energies less than 70 K, result- 
ing in a smaller output flux of cold neutrons. For 
the LANL UCN source, where substantial quan- 
tities of cold neutrons with energies around 40 
K are desired (these are subsequently cooled to 
UCN energies by inelastic scattering in solid D2), 
the absorption losses are kept acceptable by using 
4 K CH2 thicknesses in the 10-20 mm range. 

Finally, we use Figure to illustrate the dy- 
namics of neutron scattering off of moderators at 
different temperatures. Figure [7| presents the dif- 
ferential scattering cross sections for a 10 meV 
(116 K) neutron from polyethylene at 296, 77 and 
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Figure 7. Calculated differential cross section for 
scattering of a 10 meV (116 K) neutron from 
polyethylene at 296 K, 77K and 4K. 
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4 K. At 296 K, the incident neutrons with Ej = 10 
meV have energies lower than the average ther- 
mal energy of the moderator and thus most of 
the scattering amplitude is on the energy gain 
side (up-scattering). In addition, the scattering 
amplitude is almost uniform in angle; the peak at 
90 degree merely represents the cos^d^ weighting 
of the solid angle. This 47r scattering is consistent 
with the random direction of phonon propagation 
in polyethylene. As the moderator temperature 
lowers, more down-scattering occurs and the an- 
gular distribution becomes backward peaked for 
maximum energy transfer to phonons. At 77 K 
the up and down scattering are nearly equal. At 
4 K, there is almost no up-scattering and the 
down-scattering is backward peaked. This varia- 
tion in up and down-scattering is to be expected 
for neutron-moderator systems tending towards 
thermal equilibrium. 
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6. CONCLUSIONS. 

The new temperature dependent scattering 
kernel for polyethylene gives a good quantita- 
tive fit to the neutron scattering data available at 
room temperature and has a temperature depen- 
dence that is qualitatively consistent with ther- 
modynamics. The kernels are available upon re- 
quest from rhill@lanl.gov or cyliu@lanl.gov. 
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